Quarter-Filled Honeycomb lattice with a quantized Hall conductance 
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We study a generic two-dimensional hopping model on a honeycomb lattice with strong spin-orbit 
coupling, without the requirement that the half-filled lattice be a Topological Insulator. For quarter- 
(or three-quarter) filling, we show that a state with a quantized Hall conductances generically 
arises in the presence of a Zeeman field of sufficient strength. We discuss the influence of Hubbard 
interactions and argue that spontaneous ferromagnetism (which breaks time-reversal) will occur, 
leading to a quantized anomalous Hall effect. 

PACS numbers: 73.50.Jt 



Topological Insulators (TIs) are a new state of mat- 
ter with a bulk gap but protected edge/surface modes [J. 
Quantum Hall states Q are examples of time- reversal vi- 
olating TIs. Recently, time-reversal invariant TIs have 
been predicted [l[ and seen in experiments [H. They are 
characterized by a nontrivial Z 2 index in the bulk, and 
the presence of chiral edge/surface states which are ro- 
bust against localization due to static disorder. 

One of the simplest models of two-dimensional TIs is 
the Kane-Mele modelQ, which is a tight-binding model 
on a honeycomb lattice with hopping and spin-orbit inter- 
actions. A more realistic model is given by a monolayer 
of the 3D TI Alkali Iridate A 2 Ir0 3 (where A = Na, Li), 
where the various tight-binding hopping parameters are 
known from a fit to ab initio calculations 0. We will 
be working with a generic model on a honeycomb lattice 
which has no symmetries other than lattice translations, 
27r/3 rotations, and time-reversal. 

Our goal is to construct an experimentally realizable 
system with a quantized Anomalous Hall effect ( AHE) ■ 
The AHE has historically described the effect of mag- 
netic order (spontaneous or otherwise) on the Hall con- 
ductance. We will describe a proposal for observing a 
quantized Hall conductance in the absence of a perpen- 
dicular magnetic flux which can be realized, for example, 
in a monolayer of Alkali Iridate doped to |-filling, and 
subject to a Zeeman field perpendicular to the monolayer. 
Previous proposals in this direction^ have also either 
explicitly violated time-reversal, or appealed to sponta- 
neous magnetic order in interacting TIs. Our proposal 
does require strong spin-orbit couplings but does not re- 
quire the material to be a TI at half-filling. 

Generic models we consider share the following prop- 
erties: They are defined on a honeycomb lattice with two 
sites per unit cell. Including spin, there are four bands. 
In the Brillouin Zone (BZ) shown in Fig. (jTJ there are 
two kinds of special points. The T and M points are 
time-reversal invariant (modulo a reciprocal lattice vec- 
tor), while the K and K' points go into each other under 
time-reversal. Time-reversal symmetry forces the bands 
to be degenerate at the T and M a points. In a sim- 




FIG. 1: Top: The Brillouin Zone of the honeycomb lattice. 
The r and three M points are time-reversal invariant, while 
the K and K' points go into each other under time-reversal. 
The arrows show the path along which the band structure 
will be plotted. Bottom: The next-nearest neighbor hops in 
A2lrOz are shown by dot-dashed red lines that go with o\, 
the solid blue that go with with 02, and the dashed green 
that go with 173. The primitive lattice vectors are a\ and S2 
while 61, b2,bs are nearest-neighbour hopping vectors. The 
auxiliary triad shows the projection of the original crystal axes 
on to the plane of the honeycomb lattice. Note that these 
directions match the nearest neighbor hops and the labels of 
the associated a matrices. 



pie spin-independent nearest-neighbor only tight-binding 
model (e.g. Graphene) there are van Hove singularities 
at r and M a , whereas there are Dirac cones at K and 
K' . Neither of these features is generic under the addi- 
tion of time-reversal invariant spin-orbit couplings: The 
Dirac cones at the K and K' points can be gapped by 
cither a L Z S Z type spin-orbit coupling, or by a sublat- 
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tice antisymmetric potential while the van Hove sin- 
gularities are converted into Dirac cones by a Rashba 
coupling. Fig. ([2]) shows such a generic time-reversal 
invariant band structure. 

We begin with a noninteracting model with three-fold 
rotational and time-reversal symmetry, and focus on the 
generic Dirac crossings at the T and M a points. Now con- 
sider a Zceman field in the Z direction (perpendicular to 
the monolayer). This preserves the three- fold rotational 
symmetry of the lattice, but breaks time-reversal and will 
therefore provide a mass gap to the Dirac points at F and 
M a . We know that if the variation of a parameter in the 
Hamiltonian leads to the touching of two bands at a Dirac 
point for some critical value, a Chern number of ±1 must 
be exchanged between the two bands [8] as the parameter 
passes through that value Q. By the lattice rotational 
symmetry all M a must have the same exchange of Chern 
number. Since there are three M points and a single T 
point, it is clear that the Chern number exchange must 
be either ±2 or ±4. This means that upon reversing the 
Zeeman coupling, the Chern number of the lowest band 
must reverse, implying that its Chern number is cither 
±1 or ±2. Of course, the Zeeman coupling should be 
strong enough to create a hard gap (the bands should 
not overlap in energy) in order to form an insulator. We 
will show that in a simple model of A2lr0 3 with Hubbard 
interactions, such a hard gap may form spontaneously. 

To demonstrate this behaviour in detail, we start with 
the noninteracting hamiltonian for a freestanding mono- 
layer of Sodium/Lithium Iridate in the notation of Q: 
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where the spin indices have been suppressed, (ij) is a sum 
over nearest neighbors, ((ij)) is a sum over next-nearest 
neighbors, and the matrix iy in spin-space is 



U 



t' + it' a a 



(2) 



This hopping term is diagonal in the sublattice. Refer- 
ring to Fig. ([T]), the hopping is antisymmetric in the 
sublattice index (A, B). Further, each hop comes with 
a a a matrix, the index a corresponding to the projection 
of the original crystal axes on to the 111 plane. Thus, 
hopping along the primitive lattice vector ±ai carries a 
o i , hopping along the vector ±(ai — a 2 ) carries a a%, and 
hopping along ±a 2 carries a 173. Let us denote the Pauli 
matrices in the sublattice space by r a . We obtain 

h S l(k) = t' Q F (k) - t(f r (k)n - fi(k)T 2 ) - t'T 3 (F(k) ■ S) 

(3) 

where 



f(k) = f r + ifi = e lk - bl + e tk - b2 + e lk - b3 (4) 

Fo(k) = cos(fc • ai) + cos(fc • a 2 ) + cos(fc • (a 2 — ai))(5) 

F 1 (k)= sin(fc-ai) (6) 

F 2 (k)= sin(fc-(a 2 -ai)) (7) 

F 3 (k)= -sin(fc-a 2 ) (8) 



The parameters for the 111 plane of Ir in Na2lrOs 
are estimated as t = 310K, t' = -130 A' and t' = 100 K. 

One first diagonalizes F ■ a, which has eigenvalues 
±|F|, and obtains the energies at k as 



E±(k) 



t' Q F Q (k) ± sj¥\f\ 



2 + t' 2 \F\ 2 



(9) 



It is seen that the four bands come in two pairs, de- 
generate at all momenta. This is the result of inversion 
symmetry and clearly not generic. For a monolayer on a 
substrate, a Rashba spin-orbit coupling arising from the 
substrate E field perpendicular to the 111 plane is in- 
duced. A nearest-neighbor Rashba coupling leads to an 
additional term 



hn{k) = i\nG(k) ■ ctt + — i\jiG*{k) ■ <tt- 
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with cyclic permutations of 61, 62, b 3 defining G 2 , G3. 

Due to its special symmetries, this model does not 
show generic behavior for Ai?, < t'. For this parameter 
range, in addition to the usual Dirac points at T, and M a 
there are also Dirac points at K,K'. For t' > tj, (where 
t 1 * is a function of A^, but tends to 3£/8 as \r — > 0) the 
M points are unsplit. However, for Ar < t' < each 
M point splits into three Dirac points lying on the zone 
boundary. The central point is still at M and the other 
two "satellites" are symmetrically distributed about it. 
At the special point A/? = t', the satellites coalesce with 
the Dirac points at K,K'. For A# > t' the would-be 
Dirac points at K, K' are gapped, and the only Dirac 
points left are the generic ones. In the following we will 
focus on this case. 

A typical band structure including A# = OAt (satisfy- 
ing Ar > t') is shown in Fig. @. One clearly sees the 
Dirac-like crossings at the T and M a points. 

For weak Zeeman coupling, a gap will open at the 
Dirac points, but the two bands will overlap in energy (at 
different momenta) , and the system will be metallic at j 
or |-filling. Increasing the Zeeman coupling will lead to a 
hard gap and thus a quantized Hall conductance. For the 
noninteracting model, the critical Zeeman field is 0.37t. 
A calculation shows that the Chern number for these 
parameters is — 1. We note parenthetically, that even 
for the nongeneric case A/j < t' the Chern number of the 
highest energy band is nonzero (1), but different from the 
above. Once a hard gap has been established, the Chern 
number is robust against arbitrary deformations of the 
hamiltonian which do not close the gap. In particular, 
one could softly break the lattice rotational symmetry, 
add a Kane-Mele-type spin-orbit coupling or a staggered 
sublattice potential, and also tilt the Zeeman away from 
the perpendicular to the monolayer. 

Let us now consider the influence of interactions. It 
is known that a strong enough Hubard U leads to an 
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FIG. 2: The band structure for a monolayer of A^lrO^ with 
the addition of a Rashba term with \r = OAt. The path 
in the BZ followed is shown in Fig. |T]). The generic Dirac 
crossings at the V and M a points can be seen. 

antifcrromagnct in the half- filled TlQ. An on-site Hub- 
bard interaction (typically of size H- ~ 10) spontaneously 
breaks time-reversal symmetry at | filling in the Hartree- 
Fock approximation. Assuming only that lattice trans- 
lation symmetry is intact but allowing for the possibility 
of rotational and sublattice symmetry breaking, we find 
in Hartree-Fock that the magnetization is in the Z di- 
rection (perpendicular to the monolayer), and the gap 
is strongly enhanced. Because there is no spin-rotation 
symmetry in this model, the Z direction is the natu- 
ral one, and the order-parameter is Ising-like. This im- 
plies that the spin-waves are fully gapped, and that there 
should be a finite-temperature phase transition even in 
the two-dimensional material. 

Of course, it could happen that the system chooses to 
break lattice translation symmetry by forming a larger 
unit cell, as happens in Graphene at a quarter-doping jlOf . 
However, in Graphene, this is a result of Fermi-surface 
nesting, which is sensitive to the details of the band struc- 
ture, and is not generic (there are also other competing 
states in Graphene at ^-filling). Furthermore, a strong 
enough Zccman field in the Z direction will make our lat- 
tice symmetric solution energetically favorable. We will 
leave a further discussion of the issue of other possible 
ground states at quarter (or three-quarter) filling fill] to 
future work. We now speculate on the effects of long 
range (1/r) Coulomb interactions in this system, which 
can coexist with the Hubbard interaction. Long-range 
Coulomb interactions will be present because the mate- 
rial is insulating, and will likely suppress Charge Density 
Wave states, which are one avenue of translation sym- 
metry breaking. Secondly, since the band has a nonzero 
Chern number we expect skyrmionic excitations (l2j to be 



FIG. 3: The Hartree-Fock band structure at U = lOt and 
three-quarter filling in A-^IrOz with Xr = 0. It. The path 
followed in the BZ is shown in Fig. Jl}. Note the hard gap 
and the flatness of the empty band (blue). 



the lowest energy charged excitations in some regime of 
parameters. Finally, because the empty band is relatively 
flat, it may represent an environment in which analogs 
of fractional quantum Hall states will be stable when the 
band is partially occupied (l3l. Il4j. 

One can ask whether other simple lattices can sup- 
port similar states under generic conditions, where the 
only Dirac points are those mandated by time-reversal 
symmetry. The simplest square lattice two-band model 
does not support such states. The reason has to do with 
the fact that the BZ has only two M points. Repeat- 
ing the argument above Eq. ([T]), we see that the Chern 
number exchange is ±1 or ±3, implying that the states 
have Chern number ±i, zt|, which is impossible for a 2D 
band insulator. The half-filled triangular lattice with one 
band per unit cell would have a tendency to become an- 
tifcrromagnctic when interacting. Thus the honeycomb 
lattice seems to be the simplest one in which the desired 
properties can be realized. 

In summary, we have shown that in a 2D honeycomb 
lattice at ^-filling, one can obtain a band with nonzero 
Chern number with strong enough Rashba coupling, and 
an external Zceman coupling. Since the system does not 
have to be a TI at half-filling, this increases the range of 
possible material realizations. One promising material is 
A^IrO^ with A = Na, Li, with the Chern number of the 
|-filled system being ±1. In the bulk it is half-filled and 
antiferromagnetically ordered (l5j below Tjy = 15K. A 
monolayer of the 1 1 1 It plane sandwiched by two layers 
of Na would induce three-quarter doping [16[, where we 
expect it to become a ferromagnetic quantized anomalous 
Hall insulator with an Ising-likc order parameter in the 
presence of Hubbard interactions of realistic strength. It 
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may be quite difficult to grow a single monolayer doped 
in the manner described, but possibly the surface layer of 
a thin enough film will display the same properties. We 
intend to investigate this in future work. 
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